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■ ABSTRACT 

O 

' We investigate the electrodynamics of an outer gap in the meridional 



plane of the aligned-rotator. The charge depletion from the Goldrcich-Julian 
charge density causes a large electric field along the magnetic field line. The 
Q^. electrons or the positrons are accelerated by the field-aligned electric field 

^ ■ and radiate the 7-rays tangentially to the local magnetic field line. Some of 

' such 7-rays collide with X-rays to materialize as the electron-positron pairs 

^ ' on different field lines from the field line on which they were emitted. As a 

result, the electric field structure is expected to change across the field lines. 
Including these trans-field effects, we solve the formation of the electric field 
self-consistently with the curvature radiation and the pair creation processes. 
The 7-ray emission and the pair creation are treated by use of Monte Carlo 
technique. We demonstrate that the distribution of the electric field along 
the field lines is affected by both the gap geometry and the external currents 
coming into the gap through the boundaries. In the electrodynamical model, 
it has been known that the solution disappears if the current density carried 
by the electron-positron pairs produced in the gap exceeds a critical value. 
We show that the critical current density is significantly increased when the 
trans-field structure is taken into account. We also find that the location 
of the inner boundary of the gap shifts toward the stellar surface from the 
conventional null surface as the current density increases. The reason for the 
shift is derived from the stability condition of the inner boundary. We also 
argue that the ideal-MHD condition holds outside of the gap only when the 
low energy particles coexist with the high energy particles migrating from the 
gap. 
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1 INTRODUCTION 

The rotation-powered pulsars have been known as sources of pulsed radio emission. TheCompton 
Gamma-Ray Observatory (CGRO) had shown that the young pulsars are also strong 7-ray 
sources, and had detected seven 7-ray pulsars (Thompson 1999). The CGRO revealed the 
fact that 7-ray luminosity accounts for about 10% of the spin down energy. Furthermore, 
the observations for the 7-ray spectra and light curves tell us a lot about acceleration of 
particles in the pulsar magnetosphere. The energy peak in the spectra above IGeV shows 
that the electrons or the positrons are accelerated above lO^^eV. A pulsation in the observed 
hght curves imphes that the acceleration of the particles and the subsequent radiation take 
place within the light cylinder, the axial distance of which is given by = c/fi, where 
VL is the angular frequency, and c the speed of light. Although these observations provide 
constraints on proposed models, the origin of the 7-ray emission is not conclusive up to now. 

The global model for the acceleration of the particles in the magnetosphere has been 
considered as follows. The pulsar is an electric dynamo of about lO^^V. The resultant elec- 
tromagnetic energy is carried into the magnetosphere by a longitudinal current. Then, if 
there is a large electric potential drop along the magnetic field, the current carriers (prob- 
ably electrons and/or positrons) should be accelerated to high energies by the field-aligned 
electric field, and should release the electromagnetic energy by emitting 7-rays. If there were 
no electric potential drop along the magnetic field, the space charge in the magnetosphere 
should be the Goldreich-Juhan charge density (Goldreich- Julian 1969, hereafter GJ), 



where is the unit vector along the rotation axis, and B is the magnetic field. On the other 
hand, any charge depletions from the GJ value cause the field-aligned electric field E\\. 

According to the above point of view, the acceleration of particles and the subsequent 
7-ray emission have been argued with the polar cap model (Sturrock 1971; Ruderman & 
Sutherland 1975) and the outer gap model (Cheng, Ho & Ruderman 1986, here after CHR), 
the acceleration regions of which are, respectively, located near the stellar surface above the 
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magnetic poles and in the outer magnetosphere around the null charge surface {pcj = 0) 
above the last open field lines. Because the magnetic geometry of the polar caps explains 
well the observed radio polarization, it is widely accepted that the radio emission occurs 
above the polar caps. 

The outer gap model has been successful in explaining the observed light curves. By 
solving the Poisson equation for a vacuum gap, CHR concluded that the outer gap extends 
between the null surface and the light cylinder. With this vacuum gap geometry, the double 
peaks in the light curves can be interpreted as an effect of aberration and time delay of the 
emitted photons (Romani & Yadigaroglu 1995; Cheng, Ruderman & Zhang 2000). Although 
the spectrum has also been calculated with the models, in which 7-ray radiation has been 
calculated with an assumed field-aligned electric field, the models are not satisfactory in 
the sense of electrodynamics. Dyks & Rudak (2003) argued that the acceleration region 
extends to the stellar surface and the light cylinder to reproduce the observed outer-wing 
emission, and the off-pulse emission in the Crab pulsar, questioning the traditional vacuum 
gap geometry. 

A non- vacuum outer gap model was proposed by Hirotani & Shibata (1999). They have 
focused on the electrodynamics, in which screening of the electric field by the electron- 
positron pairs is taken into account. In the traditional outer gap models such as Romani & 
Yadigaroglu (1995), distribution of the electric field was not determined in a self-consistent 
manner with the radiation and pair creation processes. Hirotani & Shibata (1999) solved the 
field- aligned electric field self-consistently with the curvature radiation and the pair creation 
processes, although they worked in one-dimension along the last open field line (see also 
Hirotani 2003 for recent version of the one-dimensional model). They argued that if there 
is an injection of the particles into the gap, the position of the gap changes outwardly or 
inwardly relative to the null surface (Hirotani & Shibata 2001). 

In these one-dimensional electrodynamical models, the current carried by pairs produced 
in the gap was restricted to be about 10% of the GJ value. So, the model does not account 
for the observed 7-ray fiuxes of some pulsars, although reproducing the 7-spectra well (Hi- 
rotani, Harding & Shibata 2003; Takata, Shibata & Hirotani 2004). In the traditional outer 
gap model, on the other hand, the current densities are not calculated, but assumed to be 
the GJ values. The limitations on the current density and the 7-ray luminosity in the electro- 
dynamical model may be due to one-dimensionality. It is obvious that the electrodynamics 
is affected by the trans-field effect such that the accelerated particles causes new pairs to 
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be created on different field lines from the residing field line as pointed by CHR. Further- 
more, calculations for the light curves and the phase-resolved spectra require departure from 
the one-dimensional model. On these ground, in this paper, we extend the one-dimensional 
electro dynamical model into a two-dimensional one. 

The most important trans-field effect is caused by curved field lines. The 7-rays, which 
are emitted tangentially to the local field lines, may covert into the pairs on the different field 
lines from the field line on which they are emitted. In the present paper, to avoid additional 
complications by the three-dimension, we study a two-dimensional model. In the subsequent 
paper, we will extend the two-dimensional model to a three-dimensional one. Our aim is to 
see how the structure of the electric field, the current density, and the geometry of the gap 
are affected by the trans- field effect. 

By neglecting the magnetic component generated by the current, we assume the dipole 
magnetic field, and we consider the electrodynamics in the meridional plane of an aligned- 
rotator. This treatment simplifies the problem mathematically. Although it has been pointed 
out that an aligned-rotator is not active, that is, there are no longitudinal currents circulating 
in the magnetosphere (Michel & Li 1999), the present frame work of the aligned models will 
apply to nearly aligned-rotators with the longitudinal currents. 

We find a solution in which the particles produced in the gap carry about 30% of the GJ 
current, which is significantly larger than the value (~ 10%) obtained in the one-dimensional 
model. We also show that inner boundary of the gap shifts toward the stellar surface as the 
current increases. 

In ^21 we present the basic equations for the electrodynamics in the gap, and describe 
treatment for the pair creation process. In ^21 we show the results, and then discuss firstly the 
electric structure of vacuum and secondly non-vacuum cases. We also discuss the position of 
the inner boundary of the gap. In the final section, we compare the results with the previous 
works. We also discuss the dynamics for outside of the gap. 

In the following sections, we make statements not only about the aligned-rotators but 
also about inclined-rotators if applicable. For definite of sign of charge, we assume "parallel 
rotator" , for which the polar-caps are negatively charged. 

2 TWO-DIMENSIONAL OUTER GAP MODEL 

The present scheme consists of three main parts as follows; 
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(1) Poisson equation (§^j2^ is solved numerically with proper boundary conditions to 
give the electric field, 

(2) the continuity equation for the particles f^ ^2.2|) with pair creation term is solved on 
every field line to give the space charge density, 

(3) 7-ray field and pair creation rate are obtained by using Monte Carlo technique (§ §2.Hj) . 

These processes are iterated until all quantities are self-consistent (§ ^2.4j) . 



2.1 Poisson equation 

If stationary condition [dt + = 0, where t and are the time and the azimuth, respec- 
tively) is satisfied, that is, if there is no time variation of any quantities as seen in co-rotating 
system, then an electric field, 

E = -{ile, X r) X B/c - V$„eo. (2) 

is exerted at position r (Mestel 1999). The electric field can be expressed as the sum of the 
corotational part, —{Qe^ x r) x B/c, and the non-corotational part, — V$„co- 

Goldreich & Julian (1969) proposed that the pulsar magnetosphere is filled by plasmas, 
and that the corotational condition 

E ■ e|| = -e|| • V$„eo = (3) 

is satisfied, where e|| = B/|B| is the unit vector along the field line. If the star crust is a rigid 
rotating perfect conductor, and if the condition ^ holds all the way along the field lines 
connecting the star and the region considered, then the uniform value of $„ca over the stellar 
surface propagates into the magnetosphere. We set an arbitrary constant for $„co on the 
surface to be zero in the following. The electric field becomes E = E^o = —{^le^ x r) x B/c, 
which leads directly the corotational charge density (fT|). 

In a charge depletion region, where the condition Q does not hold, i?||-acceleration 
is indicated in the equation Q, E\\ = —e\\ ■ V$,„co 7^ 0. By substituting equation ^ in 
V ■ E = 47rp, we obtain the Poisson equation of the non-corotational potential ^nco, 

A$„co = -47r(p-pGj), (4) 

where p is the space charge density, and A is the Laplacian. 

If the variation in the azimuthal direction is negligible, we may reduce (j^ to 

Ar,9^nco = -47r(p - pcj), (5) 
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where A^^e represents (r, ^^)-parts of the Laplacian; r and 6 are the radius and the colatitude, 
respectively, in the spherical polar coordinates. This approximation is justified only if the 
gap dimension in the meridional plane is much smaller than that in the azimuthal direction. 
We must lose some effects with this simplification, but we retain the trans-field effect, on 
which we put most interest. 

We adopt an orthogonal curvilinear coordinate system based on the dipole field lines. 
The coordinates (%, C) are defined by 

T Sill ^ 

= constant along a dipole field line = x, (6) 

^Ic 

and 

T COS~^^^ 9 

= constant along a curved line (7) 

^Ic 

perpendicular to the field lines = (, 

respectively, where r is the distance from the centre of the star, and 9 is the colatitude angle 

with respect to the rotational axis. The line x = 1 draws the last open field fine, which is 

tangent to the light cylinder, and x = oo corresponds to the magnetic axis or r = oo. The 

Laplacian becomes 

_ 3 cos^ 9 + 1 4 d 3 cos^ 9 + 1 3(3cos^^ + 1) d 

'■'^ ~ ^ Xsm^9d^ ^ 4cos3^ ^ 4Ccos3^ ^C' ^ ' 

The non-corotational electric field E„co is obtained by 

^ d^nco ^ V3 COS^ 9 + 1 d<^nco 

II 2zui,cosy^9 dC ' ^ ^ 

os± zuic sm 9 ox 

where ds\\ — 2 cos^/^ 9/\/2,cos^ 9 + IdC, and ds±_ — s\v? 9/\/Zcos^ 9 + Idx are the line ele- 
ments along the field lines and the perpendicular curved lines, respectively. Since dipole 
field is curl-free, we use pcj ~ —QBz/2nc as the GJ charge density. 



2.2 Continuity equations for particles 

Let us denote velocity v as 

v = -yiieii +Vd, (11) 

where is the drift velocity owing to the electric field in the gap. The drift velocity is 
separated into the corotational and the non-corotational parts. 
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Vd = Vco + C — , (12) 

where v^o = ^^e^ x r. In the later section O, we shall see that the magnitude of the 
non-corotational drift is negligibly small in comparison with the corotational drift. By using 
V ■ (nv) = B ■ V{v\\n/B) + Vldn/dcj) and imposing the stationary condition [dt + VLd^f, = 0), 
the continuity equations become 

B.v(^)=±S(r). (13) 

where S{r) is the source term due to the pair creation, and denote the number density 
of outwardly and inwardly moving particles, respectively. We suppose that the accelerated 
particles have almost the speed of light, ~ Vu+v'^^ ~ c^, such that the parallel component 



is assumed to be f|| = yc^ — {-cu^l)'^. With this parallel motion, the continuity equations 
(|T3|) yield the particle's number density, provided that the source term S{r) is given. 

If one assumes that the particle's motion immediately saturates in balance between the 
electric and the radiation reaction forces, the Lorentz factor of the particles within the gap 
is simply given by 

Wr) = + ij , (14) 

where Rc is the curvature radius of the magnetic field line. This saturation simplifies the 
problem significantly. However, the particles in the gap will have the saturated Lorentz factor 
Tsat only when typical accelerating time tac, radiation damping time t^ and crossing time 
tcr = W\\/c, where W\\ is the gap width along the magnetic field lines, satisfy the condition, 

tcr > tac and tcr > td- (15) 

Hirotani, Harding & Shibata (2003), who solved unsaturated motion of the particles, showed 
that the condition tcr ^ tac is not satisfied effectively for some pulsars. Although the un- 
saturated motion affects the 7-ray spectrum below IGeV (Hirotani et al. 2003; Takata et 
al. 2004), effect on the electric field is less important. As far as the electrodynamics are 
concerned, the assumption of the saturation will provide a good approximation. 
The source term can be expressed as 

S{r)= de^[rip+G+ + r]p_G^], (16) 

where is the photon energy normalized by meC^, (7+ and denote the distribution 
function of outwardly and inwardly propagating 7-rays. The pair creation rate rjp per unit 
time per 7-ray photon with the energy of is given by 
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Figure 1. The pair creation cascade model in tlie meridional plane. The 7-rays are radiated by the curvature process, and are 
beamed in the direction of local magnetic field. The radiated 7-rays may convert into the pairs by the pair creation process. The 
Xl) Xupi Cini and (^out present the lower, upper, inner, and outer boundaries of the gap, respectively. 



. , dNx 



(17) 



Hh dex 

where dex ■ dNx /dex is the X-ray number density between energies meC^ex and meC^{ex + 
dex), cos~^ /ic is the colhsion angle between a X-ray and a 7-ray, rrieC^eth = 2meC^ /{I — fic)e^ 
is the threshold X-ray energy for the pair creation, and cxp is the pair creation cross-section, 
which is given by 
3 



where 

v{e^,ex) 



16 



(Tyfl — V 



(18) 



'1 



1 



/ic e^ex 



and (Jt is the Thomson cross section. 

In this paper, we take the blackbody radiation from the star as the X-ray field. At the 
distance r, the photon number density is given by the Planck law. 



dNx 1 f2TTmeC^Y f \ 



X 



^ (19) 

dex 4:71 y ch J KAirr'^ J exp^rriec'^ex / kTs) — l' 

where As is the emitting area, and kTs refers to the surface temperature. For the values of 
As and T^, observed ones may be used. 



2.3 Pair creation model 



We calculate ^(r) by simulating the pair creation process with Monte Carlo method. 

FigH shows the pair creation cascade model schematically. Accelerated particles emit 
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7-rays tangentially to the local field lines by the curvature process. The 7-rays deviate 
from the field line on which they are emitted. A small fraction of the emitted 7-rays may 
collide with the X-ray photons to materialize as electron-positron pairs. For example, Go 
photons emitted at position Xq in FigUl materialize as SNi, 6N2 ■ ■ ■ Go) pairs at positions 
Xi,X2, ■ ■ ■, respectively. The new born electron-positron pairs such as the pairs depicted at 
the position Xi and X2, will be separated immediately by the field-aligned electric field 
after their birth. The separated pairs also radiate the 7-rays which become new seeds of a 
sequence of pair creations. The separated pairs carry a current, and provide a space charge 
to screen out the electric field. 

On the other hand, as would be expected near the upper boundary of the gap, if the 
potential drop along the field lines is too small to separate a pair, then the pair would go 
out the gap together, such as pair depicted at position X3. This kind of pairs does not affect 
the electrodynamics in the gap, in spite of the fact that they are produced in the gap. The 
pairs produced outside of the gap do not affect the dynamics either, because they do not 
directly come into the gap. 

To obtain the source term S{r) at each position, we calculate the annihilated photon 
number 6N at the position r. In terms of mean free path of a 7-ray with energy e^, lp{r, e^) = 
c/rip{r,e^), the number of annihilated 7-rays (produced pairs) between the distances s and 
s + ds from the emission point (s = 0) is 



For the typical pulsar parameters, the mean free path Ip is found to be much longer than 
the light radius, Ip ^ wic, so that we can approximate equation (pUj) as 6N ~ {Go/lp)ds. 
We assume that Ip does not depend on the position, taking the value at r = 0.5wic with the 
collision angle 7r/2 {fic = 0). This assumption is probably a rough treatment because the 
X-ray number density ()19p decreases with the distance from the centre of the star. In the 
present model, the pair creation position is given by uniform random number for the given 
element of the 7-ray flux. It is notable that the 7-rays are beamed in the direction of local 
particle's motion v = V|| -|- Vco^^, so that an angle between the direction of the propagation 



and the meridional plane is 6^ = cos^^[y1 — {^Iw/wicy]. Hence, the displacement in the 

meridional plane is reduced to be scos6'-y, where s is the three-dimensional path length of 
the 7-rays. 

Using the above pair creation model, the source term S{r) is obtained by the following 




(20) 
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procedure (see also FigH]). 1) We calculate the number Gq of the emitted photons per unit 
time per unit volume at the position Xq. 2) We calculate the number of photons converted 
into pairs in a distance so(~ voic) by 5G = Go[l — exp{— Jq° 1/lpds')]. 3) We divide the 6G 
photons into n flux elements, and pick up n points on the three-dimensional path length 
{si < So,i = l,---n) randomly. The pair creation point Xj in the meridional plane are 
deduced from |Xj — Xo| = SjCos6'^. 4) We calculate the produced pairs per unit time per 
unit volume as the source term ^(Xj). If potential drop between a pair creation point and 
boundaries is small (e5$„co < lO'^mgC^ in this paper), we do not count the pair as the source. 
The emitting points of the 7-rays are scanned along the field line (from the inner boundary 
to the outer boundary) and across the magnetic field lines (from the lower boundary to the 
upper boundary). This direction of scanning across the field lines has an advantage because 
the 7-rays propagate always in convex side of the field lines. The particle's number density 
along the magnetic field line follows the equation (|T3|) . For the simplest treatment, all the 
curvature photons are assumed to be emitted at the critical energy, 

. ±^ . 0.1 (^) GeV, (21) 



477 Rc ■ VlOOs-V VlOV 
with the rate, 

9 hRc VlOV \wiJ ^ ' 

We introduce the conventional dimensionless variables; 



r=— r, ^p = W 7^ , 23) 

$nco = -^Kco, E = -V<lnco, (24) 

and 

(3 = -, (26) 

c 

where Bic is the strength of the magnetic field at the point on which the last open field line 
is tangent to the light cylinder. By using dimensionless charge density, 
P ^ 

P=mj2^c = P^,^'^-'~^^ PC. = -5., (27) 
the Poisson equation (@)) and the continuity equation (|13|) are rewritten as 

A$„eo = -[|-(j+-J-) + 5,], (28) 

PI I 



A pulsar outer gap model with trans-field structure 11 



and 



B ■ Vj± = ±^(f) 



(29) 



respectively. 

2.4 Boundary conditions and method of calculation 

We assume that the lower and the upper boundaries are laid on the magnetic surfaces labeled 
by xi and Xup, respectively. The inner and the outer boundaries are defined by the surface 
on which E\\ vanishes, 



where Cm and (out, which are in general function of denote the location of the inner and 
the outer boundaries, respectively (Fig|T)). 

The conventional polar cap and outer gap models have their current densities flowing 
in opposite directions to each other. For this difference, it is impossible that these two 
acceleration regions are located on the same magnetic fleld lines. Therefore, we anticipate 
that the inner, the upper and the lower boundaries are directly linked with the star without 
any potential drop. We then impose 



The conditions (jHUj) and (jHT|) are not satisfled on arbitrary given boundaries, because both the 
Dirichlet-type and the Neumann-type conditions are imposed on the inner boundary Qnix)- 
We, therefore, solve the location of the inner boundary. By moving the inner boundary step 
by step iteratively, we seek for the boundary that satisfies the required conditions. 

We have to consider the stability of the gap. The stability condition for the inner bound- 
ary is dE\\/ds\\ ^ such that external electrons may not freely come into the gap through 
the inner boundary. This is similar to the force-free surface in the quiet aligned model by 
Michel (1979). In the same way, for the outer boundary it is dE\\/ds\\ ^ 0. Furthermore, 
we postulate that the gap is unstable if the field-aligned electric field changes its direction. 
In other word, the sign of the field-aligned electric field is positive definite, as it guarantee 
that the positrons and the electrons propagate outwardly and inwardly, respectively. For the 
non-vacuum case, if super-GJ (|p| > |pgj|) region is dominant in the gap, then the electric 
field may change its direction. 



E||(c.„) = ^ii(Co«t) = 0, 



(30) 



'^nco(Cin) '^nco(Xup) '^nco(Xi) 0. 



(31) 
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The continuity equation ()29|) satisfies conservation law of the longitudinal current density, 
that is, 

itot = j+ + 3~ = const along field line, (32) 

where jtot is the total current density in units of the GJ value. The total current is made up 
of the three components, 

(1) jg, the current carried by the electrons and the positrons produced in the gap, 

(2) ji, the current carried by the positrons (e.g. originating in sparking on the stellar 
surface) coming into the gap through the inner boundary, 

(3) j2, the current carried by the electrons (e.g. originating in pulsar wind region) coming 
into the gap through the outer boundary. 

We cannot deduce the value of jtot from any gap model, because the current circulates in 
the magnetosphere globally. The outer gap, the polar cap and the pulsar wind may interact 
each other through the current, so that the current should be determined by some global 
conditions. In our local model, therefore, we regard the values of (jg, ji, ^2) as a set of free 
parameters. The boundary conditions on the current at the inner and the outer boundaries 
are given by 

j+iCin) = jiix), (33) 
and 

J-{Cout)=j2{x), (34) 
respectively. In terms of {jg,ji,j2), the conservation low (fH^ yields 

ji + 32 + jg = jtot- (35) 

As has been mentioned, any local model should have some degrees of freedom because 
the model includes global quantities such as the electric current. The free parameters cannot 
be eliminated unless one solves the interaction with other parts of the magnetosphere (e.g. 
with pulsar wind). Furthermore, we have the free boundary problem with both the Dirichlet- 
type and the Neumann-type conditions. Hence, choice of a set of free model parameters is 
somewhat arbitrary. For convenience for the iteration, we choose the model parameters as 
follows. The external current sources (ji, j2) are obviously suitable to be free (e.g., in the 
traditional gap models, one assumes ji = j2 = 0). The lower boundary would be chosen as a 
model parameter, but it is traditionally set to be the last open field line, that is, xi = 1- We 
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give also the locations of the upper boundary Xup and the outer boundary (out{x)y thereby 
we obtain a specific outer gap, namely we solve jg{x) and the location of the inner boundary 
Cinix)- Giving different (outix) and Xup, we obtain a series of models for which jg{x) and 
Cm(x) are obtained. 

One specific solution is determined by iteratively as follows. 1) By assuming the field- 
aligned electric field and the position of inner boundary Cm(x)) we start our calculation. 2) 
We calculate the particle's Lorentz factor (jl4|) from E^^. 3) We solve the curvature radiation 
and the pair creation processes by the method described in ^ ^2.31 and then obtain the 
particle's number density in the gap. 4) Having the space charge density, we solve the Poisson 
equation to obtain new field-aligned electric field and the inner boundary that satisfies the 
required conditions. 5) We repeat 2)-4) until all the physical quantities are self-consistent. 

The values of all the physical quantities go quickly toward convergence through the it- 
eration. For example, if the field-aligned electric field is over estimated, then the resultant 
particle's number (or charge) density calculated with the radiation and the pair creation pro- 
cesses becomes larger than the solution. With this over estimated number density, however, 
the Poisson equation gives a smaller electric field in the next step. 

3 RESULTS 

3.1 Model parameters 

As described in the last section, we take (ji, j^), {Xh Xup) and (out as the model parameters. 
We fix X; = 1 (i.e. the lower boundary is the last open field line), and (out = 1-6 such that 
the trans-field path length of 7-rays in the gap may becomes significant. To examine the 
effect of variation in the external current, we consider the following four cases; 

Case 1: (ji,j2) = (0.0,0.0) for Xup ^ X > x', but (ji,j2) = (10-^ 10-^) for x' > X ^ Xi, 
where x' = Xi + 0.02(x«p - Xi), 

Case 2: (ji, ^2) = (0.1,0.1) for Xup > X > Xh 
Case 3: (ji, ^2) = (0.2,0.0) for Xup>X> Xh 
Case 4: {3^,32) = (0.0,0.2) for Xup>X> Xi- 

In Case 2, the positronic and the electronic currents of 10% of the typical GJ value are 
injected into the gap from the inner and the outer boundaries, respectively. In Case 3 (or 
Case 4), on the other hand, the only positronic (or electronic) current comes into the gap 
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through either side of the gap. In Case 1, we consider the case of no external current, that is, 
the whole current is carried by the pairs produced in the gap. Because we need seed 7-rays 
to start the pair creation cascade in Case 1, we assume a very little currents (ji = j2 = 10~^) 
in a thin layer (2% of the trans-field thickness) just above the last open field line. 

In general, jg increases with Xup (i-e. with the trans-field thickness). But, if jg becomes 
larger than a critical value, we cannot get the stable solutions because the field-aligned 
electric field changes its sign in the gap. In this paper, we regard the maximized Xup as 
a natural boundary, and thereby we obtain Xup ~ 1-22 (Case 1), ~ 1.2 (Case 2), ~ 1.13 
(Case 3) and ~ 1.21 (Case 4). For the pulsar parameters, we have used the values of the 
Vela pulsar; Vt = 70.6s-\ B, = 6.68 • IO^^q^ j^j^ = 150eV, and As/A, = 0.066 (Ogelman, 
Finley & Zimmerman 1993), where Bg is the polar field strength, and A^, is the area of the 
stellar surface. 

3.2 Vacuum gap 

Before we investigate the electrodynamics of the non-vacuum gap, let us see some properties 
of the vacuum solution. 

If the gap width along the magnetic field lines is much shorter than the trans- 
field thickness D_i_, then the one- dimensional approximation works out well. From the one- 
dimensional Poisson equation, dE\\/ds\\ = B^, the inner and the outer boundaries satisfy 



where B^ is the dimensionless GJ charge density, s\\in and s\\out are the locations of the inner 
and the outer boundaries, respectively, in terms of the arc length along the field line from 
the stellar surface. In this case, the field-aligned electric field E\\ varies quadratically along 
the field hues such as one displayed in FigO 

CHR pointed out that the trans-field thickness is suppressed by the pair creation, and 
postulated that the gap width is much longer the trans-field thickness, W\\ ^ D±. By solving 
the Poisson equation with this slab-like geometry, they found that the inner boundary is 
located close to the null surface. While CHR solve the Poisson equation with a recti-linear 
geometry, the present model gives a vacuum solution with the dipole magnetic field in 
Fig|2fa). The result retains a property of CHR; that is, the inner boundary is located near 
the null surface. This property does not depend on the location of the outer boundary. 




(36) 
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Figure 2. A solution of the vacuum gap. (a) The gap in the magnetosphere (filled region). The solid and dashed lines show 
the dipole field lines and the curve lines that perpendicular to the field lines, respectively. The dotted line is the null surface, 
(b) The solid, dashed, and dotted lines show the distributions of the GJ density, the two-dimensional term, and acceleration 
field, respectively, along a magnetic field line. 



In the elongated gap case, instead of the equation (jHSl), the inner and the outer boundaries 
satisfy 



0, 



(37) 



where Aj^ is the trans- field part of the Laplacian. Figl^b) shows the distribution of B^, the 
transverse term A±^nco, and E\\ on a magnetic field. The transverse term represents the 
effect of surface charge on the upper and the lower boundaries. From FigEl^b), we find that 
the GJ charge density and the transverse term almost balance out in most part of the gap, 
Bz ~ — Ax$„co- Therefore, the distribution of E\\ has a plateau in the middle of the gap as 
the dotted line in FigI2fb) shows. 



3.3 Non- vacuum gap 

In this subsection, we discuss the electrodynamics of the non-vacuum gap. The lowest order of 
the Fourier components of the two-dimensional term is represented by A±^nco ~ —^nco/D±, 
and therefore the two-dimensional term plays the role of a positive charge density in the 
gap. 



3. 3. 1 Electric field structure 

Fig El shows the distributions of the field- aligned electric field E\\ along the magnetic field 
lines. The solid, dotted and dashed lines represent, respectively, the electric field on the 
magnetic field lines threading the gap at 25, 50, and 75 % of the trans-field thickness from 
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the lower boundary. The abscissa refers to the arc length from the null surface. The positive 
and the negative values indicate outside and inside, respectively, with respect to the null 
surface. 

The variations of E\\ along the magnetic field line for all cases are similar to that for the 
vacuum case displayed in FigEl In these cases, because the trans-field thickness is smaller 
than the gap width, the transverse term Ax$„co is important. The value of E\\ is determined 
in such a way that Aj_$„co; the local charge density p and the GJ charge density —B^ 
almost balance out in most part of the gap, —K^^nco — P ^ (see FigEj). As a result, the 
distribution of the electric field shows a plateau in contrast with the quadratic shape (such 
as FigCj) in the one-dimensional model, in which the transverse term is ignored. 

The effect of the external currents (ji, J2) produces the difference in the distribution of 
E\\. For example, E\\ on the plateau in Case 3 gradually decreases with the distance, while 
that in Case 4 gradually increases. In terms of ji and 72, the local charge density ipTj) is 
rewritten as p = B{ji — j2 — jg + 25j)//5||, where 5j represents the current density carried by 
the positrons produced between the inner boundary and the position considered (e.g. 5] = 
at the inner boundary, and 6j = jg at the outer boundary). For Case 3 (ji 7^ 0,^2 = 0), the 
charge density at the outer boundary is 13{ji+jg)/ indicating that the external positronic 
current ji assists the current jg in screening of E\\ at the outer boundary. Therefore, the 
strength of E\\ near the outer boundary is smaller than that near the inner boundary. In 
the same way, the external electric current j2 for Case 4 assists the screening of E\\ at 
the inner boundary, where the charge density is -B(— J2 ~ jg)/ Hence the strength of E\\ 
near the inner boundary is larger than that near the outer boundary. To be expected, if 
both the positronic and the electronic currents are injected the respective boundaries such 
as Case 2, these two contributions in the Poisson equation cancel each other, so that the 
plateau becomes almost flat. 

Another effect of the external currents (ji, j2) appears in the position of the inner bound- 
ary. The inner boundary for Case 3 is located outside of the null surface, and as a result 
the gap does not include the null surface, although the boundaries for the other cases are 
located inside of the null surface. This difference is discussed in ^ ^3.3.31 
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3.3.2 Current density 

The trans-field structures found in the current density and the field-aligned electric field E\\ 
are summarized in Fig|3 The solid lines show the distribution of the current jg across the 
field lines. The maximum value (arbitrary units) of E\\ on each field line is also plotted in 
FigEl 

We find that maximum value of jg for each case is jmax ~ 0.2 (Case 1), ~ 0.17 (Case 2), 
~ 0.07 (Case 3) and ~ 0.3 (Case 4). The total current density on each magnetic field line is 
represented by jtot = jg + ji + j2- In Case 1, for example, because no particles come into the 
gap through both the boundaries, (i.e., ji = j2 = 0), the total current is carried only by the 
pairs produced in the gap, that is, jtot = jg- The total current densities at the maximum are 

Jtot Jmax 

~ 0.2 (Case 1), j^ax + ji + j2 ~ 0.37 (Case 2), jmax + ji ~ 0.17 (Case 3), and 
jmax + j2 ~ 0.4 (Case 4). 

As mentioned in ^ ^3.11 the solution disappears if jg exceeds the critical value because the 
super-GJ region becomes dominant in the gap so that the field-aligned electric field changes 
its direction. However, it is noteworthy that the values jg for Case 1, Case 2 and especially 
Case 4 exceed the critical value (~ 0.1) of the previous one-dimensional model. As an effect 
of the trans-field structure, the transverse term A^$„co behaves as the positive charge. 
Therefore, if the super-GJ points appear in the negatively charged region, the effect of it 
in the Poisson equation can be reduced by the transverse term. In Case 4, for example, the 
super-GJ points tends to appear in the negatively charged region around the inner boundary 
due to the external electric current j2 (see Fig^. Then, we find that the shape of its inner 
boundary (see next subsection) and the gap structure are formed so that the transverse 
term cancels the super-GJ points as Fig0] shows. Because we increase jg until the transverse 
term marginally sustain the super-GJ region, we obtain the critical value (~ 0.3) of Case 4 
increasing significantly from the one-dimensional one (~ 0.1). The situation such as Case 4 
is very important for obtaining the solution of the increased current. On the contrary, if the 
super-GJ space charge density appears near the outer boundary, it is positive, so that the 
transverse term does not reduce the super-GJ charge density but the opposite. In such a 
case (e.g. Case 3 with positive external current), the critical current density is rather small. 

In the outer gap model by CHR, the pairs are supposed to be created much in the 
convex side (i.e. in the upper half), and a large current fiows there, while they assumed 
the vacuum region in the lower half. This type of gap is realized in Case 1 by assuming no 
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external currents. From the panel for Case 1 in Figl^l we find that the vacuum-like region 
occupies a lower half of the gap, and the current flowing (non- vacuum) region is restricted in 
the upper half region. As far as the trans-field structure is concerned, our electrodynamical 
model reproduces the structure similar to that assumed in CHR. However, the location of the 
inner boundary manifest itself in different manner, as shall be shown in the next subsection. 

FiglHl shows that the current jg increases with increasing height, and then decreases. The 
increase is simply because the 7-rays propagate the convex side of the field lines and make 
pairs. The decrease is caused by the following effects. Firstly, the upper part is illuminated 
mainly by the 7-rays which were radiated in the middle region, while almost of the 7-rays 
emitted at the lower half part escape through the inner or the outer boundaries before 
reaching the upper region. Secondly, the field-aligned electric fields, the Lorentz factor (oc 
E]^'!^) and in turn the energy of the emitted 7-rays (oc E^'l^) decrease with increasing of the 
height in the upper half region. Furthermore, the threshold energy of the X-rays colliding 
with the low-energy 7-rays moves to the Wien region of the Planck law, and thereby, the 
mean free path of the 7-rays sharply increase with decreasing of its energy. By combinations 
of these effects, the current density jg decreased with increasing of the height in the upper 
half region. 

In contrast, the exponential increase of the currents near the upper boundary are pre- 
dicted in CHR. This case may be obtained, if one takes non-thermal X-ray fields with a 
modest power-low index such as the Crab pulsar, or if one considers a case in which the 
strength of E\\ is ineffective to the Lorentz factor and the energy of the 7-rays (e.g. the un- 
saturated motion of particles). These possibilities are investigated in the subsequent papers. 



3. 3. 3 Location of the inner boundary 

The filled regions in FiglHl show the shape of the gap. The heak-Wke shape of the inner 
boundary is remarkable when one compares it with that of the vacuum case. The magnetic 
field line on which the current density jg is maximized at jmax runs through the cusp of the 
beak. 

An important feature of the external currents (ji, J2) appears in the position of the inner 
boundary of the lower end part, where few pairs are produced. If ji and j2 are equal such as 
Case 1 and Case 2, the inner boundary is located close to the null surface. If ji and j2 are 
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Figure 3. The electric field structure for Case 1, Case 2, Case 3, and Case 4. The solid, dotted, and dashed lines show the 
distributions of the field-aligned electric field on the three different magnetic field lines, which locate, respectively, 25, 50, and 
75% of the gap trans-field thickness measured from the last open field line. The abscissa represent the arc length with the origin 
at the null surface. 
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Figure 4. The distributions of the GJ charge density (solid line), the transverse term (dashed line) and the space charge 
density (dotted line) Oil the nia,gii6tic fislci line, on wliicli jg jmax 

for Case 4. 



not equal such as Case 3 and Case 4, on the other hand, the inner boundary shifts outside 
or inside. This is because the effective charge neutral shifts relative to the conventional null 
surface by the space charge caused by the external current. This effect is originally pointed 

out by the onc-dimcnsional model of Hirotaiii & Shibata (2001). 

As mentioned above, the inner boundary at the lower part region of Case 1 comes to 
close the null surface. On the other hand, we find that the inner boundary at the cusp. 
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distribution of the current jg carried by the pairs produced in the gap. The dashed line shows the maximum value (arbitrary 
units) of on each field lines. 
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Figure 6. The geometry of the outer gap (filled region in each panel). The solid and dashed, and dotted lines are the same as 
Fig'El On the dashed-dotted line in each panel, the condition, jmax +32 — ji = B^/B, is satisfied (see text). 



where jg = jmax, leaves from the null surface and moves toward the stellar surface. Let us 
investigate this displacement. From the stability condition discussed in § ^2.4| we have 

dEu 



ds\ 



Ax$„co + [B{u - j-)/P\\ + B,] ^ 0. 



With the help of BJB ^ 1 and the fact that both A_l$, 
the cusp, we have 



(38) 

at the and ~ 1 hold on 
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Umax + J2 - Jl) ^ 5,(Ci„)/5(Ci„) ^ 1, (39) 

where we have used the fact that j+ — j„ becomes ji — j2 — jmax at the cusp. This condition 
indicates that the inner boundary is located between the stellar surface where B^/B = 1 and 
the position on which the local charge density caused by the longitudinal currents equals 
the GJ charge density. 

In Case 1 and Case 2, because ji equals j2, the condition (jHHj) becomes 

jmax ^ B,{U)IB{U) ^ 1- (40) 

The dashed-dotted lines for Case 1 and Case 2 in FiglHl show the positions on which the 
condition jmax = Bz/B {jmax = 0.2 in Case 1 and jmax = 0.17 in Case 2) is satisfied. We find 
that the cusp is located significantly closer to the position where the space charge density, 
—jmaxB, equals the GJ charge density, —B^. Noting that B^/B increases with decreasing 
of r, we expect that the cusp shifts toward stellar surface as jmax increases. If the current 
carried by the particles produced in the gap were the GJ value (i.e. if jg = 1), the cusp 
would be located on the stellar surface, and therefore the outer gap would start from the 
stellar surface. However, because we have no solution if jg exceeds a critical value (~ 0.2 
for Case 1), we do not obtain the solution in which the inner boundary touches the stellar 
surface. 

For Case 3 and Case 4, the condition (j39|) is rewritten as (jmax—ji) ^ Bz{Cin)/B{Qn) ^ 1 
and {jmax +72) ^ Bz{(in)/ B{(in) ^ 1, respectively. The dashed-dotted lines in the lower 
panels of FiglHlshow the positions which satisfy the condition {jmax — ji)B = B^ for Case 3 
and {jmax + j2)B = Bz for Case 4. Such as Case 1 and Case 2, the cusps in these cases also 
take the position on which the space charge density equals the GJ charge density. Thus, we 
find that the outer gap extends from the point on which the charge density caused by the 
longitudinal current made up of {jmax, ji, j2) equals the GJ charge density. 



4 DISCUSSION 

We have studied the electrodynamics of the outer gap including the trans-field structure 
in the meridional surface. We have found that the inner boundary shifts toward the stellar 
surface from the null surface as jg or j2 increases, and forms the beak shape, which was 
found in the vacuum gap geometry. We have shown the specific solution, of which the 
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current density jg carried by the particles produced in the gap becomes significantly lager 
than the upper limit given by the one-dimensional model {jg ~ 0.1). 

4.1 Comparison with previous works 

After CHR, the traditional shape of the outer gap have been based on the vacuum model, 
in which the gap extends between the null surface and the light cylinder. Although the tra- 
ditional outer gap models have been successful in explaining overall features of the observed 
light curves (e.g. double-peaks in one period, and the presence of the bridge emission), the 
model cannot account for the presence of the outer-wing emission and of the off-pulse emis- 
sion from the Crab pulsar. Recently, the caustic model by Dyks & Rudak (2003), and Dyks, 
Harding & Rudak (2004) reproduces the observed light curves better than the traditional 
outer gap model. The caustic model demanded that the emission region extends from the 
polar caps to the light cylinder. But, the reason that the gap extends on both sides of the 
null surface is not evident. Our two-dimensional electrodynamical model indicates that if 
a marginal stable condition, {jg + j2 — ji)-B(Cm) = Bz{Cin) holds, then the inner boundary 
shifts toward the stellar surface as jg and j2 increase. 

We have actually demonstrated that the gap extends on both sides of null surface in 
Case 4 (FigjHl), in which the electrons come into the gap through the outer boundary. In the 
global magnetospheric models (e.g. Mestel et al. 1985, Shibata 1995), the electrons circulate 
the magnetosphere globally and turn back into the outer gap. Although the observed light 
curves may be reproduced better by assuming the injection of the electrons at the outer 
boundary (ji = 0,^2 7^ 0), it is not certain whether the injection from the outer boundary is 
favored for reproducing the 7-ray spectrum. This will be studied in a subsequent paper. It 
is notable that, in the one-dimensional model, the observed spectra have been reproduced 
by assuming the injected positrons from the inner boundary (Hirotani et al. 2003; Takata 
et al. 2004). 

Although the present model is very simple, we must argue that two or three dimensional 
model will open a way to diagnose the current system in the outer gap when the observed 
phase resolved 7-ray spectra and light curves are compared with the model predictions. The 
external currents (ji,j2) strongly control the 7-ray spectrum by changing the electric field 
strength. On the other hand, the current density jg and its distribution across the field lines 
reflect the pulse shapes through the location of the gap boundary. Thus, combining the 
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spectra and the pulse shapes, the electrodynamical model enables us to know the current 
system in the outer gap. 



4.2 Outside of the gap 

In this subsection, we discuss the dynamics in outside of the gap. We shall show that co- 
existence of low energy particles and the high energy particles streaming from the gap is 
required to hold the ideal-MHD condition, E\\ = 0, in the outside of the gap. We shall lead 
this conclusion with the reductio ad absurdum. 

Taking Case 1 (ji = j2 = 0), let us consider the region between the stellar surface and 
the inner boundary. We assume at the moment that only the electrons escaping from the gap 
exist in that region and migrate toward the stellar surface. The Poisson equation and the 
energy conservation for the electrons in the one- dimensional form are 

^ = -I^-M + 5 (41) 

and 

|^ = -^||, (42) 
ds\\ 

respectively, where S|| is the arc length along the field line from the stellar surface, Tq is the 



Lorentz factor of the particles on the inner boundary, and F = 1/ y^l — (3'^^ — The first 
term in the right hand side of the equation (j^T|l represents the two-dimensional effect. 

Near the inner boundary but outside of the gap, because ~ 1 and because F ~ Fq, the 
equation (PT|) can be rewritten as dE\\/ds\\ = —jgB + Bz- Since jgB ^ i?^ is satisfied, we find 
dE\\/ds\\ ^ 0, where the equality is satisfied on the inner boundary. By noting E\\ = on the 
boundary by the definition, the field-aligned electric field has the negative sign in the outside 
of the gap, so that the Lorentz factor of the escaping electrons decreases. If the decrease of F 
from Fo becomes a considerable value, the transverse term dominates the right hand side of 
equation ()41|) . This transverse term causes further decrease of and subsequent exponential 
development of the Lorentz factor. Thus, we find that the ideal-MHD condition does not 
hold between the inner boundary and the stellar surface. On these ground, we conclude that 
other electrons populate between the stellar surface and the inner boundary to hold the 
condition E\\ = (i.e. the ideal-MHD condition). These electrons will be supplied from the 
stellar surface or the electronic cloud on the polar caps by the negative field-aligned electric 
field, which will appear unless there are such electrons. 
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The same argument is applied for outside of the outer boundary, and therefore the ions 
(or positrons) in addition to the positrons escaping from the gap are required. 

4.3 Displacement of the outer boundary 

In all the four cases in §21 we assume that each position of the outer boundary is located 
far outside from the null surface. Then expansion of the gap in the trans-field direction is 
suppressed by the pair creation. As a result, the trans- field thickness becomes shorter 
than the gap width W\\. 

In our model, the location of the outer boundary Cout is free (^ §2.4|1 . but determine the 
value of the total current running through the gap. We may put the outer boundary near the 
null surface. In such cases, because W\\ is short, we have a gap which has a larger than 
those of Case l~Case 4. In this wedge-Wke gap, because the trans-field propagating distance 
of the 7-rays in the gap is negligible as compared with D^, the magnetic field lines are 
effectively straight in the gap. We then anticipate that the solution of our two-dimensional 
model with this geometry is effectively the same as that of the one-dimensional model. 

Fig 13 shows a solution of the wedge-hke gap, where the model parameters ji =32 = 0.1, 
which are the same as Case 2. As a result, there is no solution if jg becomes greater than 0.1, 
which is smaller than the critical value of ~ 0.17 for Case 2. Furthermore, the distribution 
of the field-aligned electric field is quadratically, and is difference from one of the slab-like 
gap such as Case l~Case 4. 

4.4 Drift motion 

Let us examine validity of the assumption that the non-corotational drift motion is negligible 
as compared with the corotational motion. From equation (fT^ . the drift motion around the 
rotational axis may be expressed as 

= zuQ + c— , (43) 

where E± = —d^nco/ds±. The first and the second terms of the right hand side of the 
equation ()43p represent the corotational and the non-corotational motion, respectively. Fig|Hl 
shows the drift velocity of the non-corotational part evaluated on the outer boundary. The 
plus (or minus) sign for the abscissa represents the direction (or counter-direction) of the 
stellar rotation. From this figure, we find that the drift motion becomes super-corotation in 
the lower part of the gap, and sub-corotation in the upper part. In our calculation, we find 
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Figure 7. A solution of tlie wedge-like gap. The external currents (j'l, ^2) = (0.1, 0.2) is the same as Case 2. (a) The geometry 
of the outer gap (dark region). The solid, dashed, and dotted lines are the same as in Fig|E| (b) The distribution of the field- 
aligned electric field. The solid, dotted, and dashed lines are the same as in Fig|21 (c) The trans field structures of the current 
and the field-aligned electric field. The solid and dashed lines are the same as in FiglSl 



that non-corotational motion is a few per cent of the corotational motion (i.e. |Eco| ^ 
The rotational period P3 of drift motion of plasmas due to E\\ in the outer gap is suggested 
to be P3 ~ 13Pi, where Pi is the corotational period, and we used the maximum drift 
velocity |f | ~ 0.06c at zu = O.Szuic- Since we have neglected the structure of the azimuthal 
direction, we cannot calculate the drift velocity due to the azimuthal electric field. However, 
we may anticipate that the strength of the azimuthal electric field are same order (or less) 
as compared with the values of the trans- field electric field E±, unless the dimension in the 
azimuthal direction is much smaller than that in the meridional plane. Therefore, we safely 
neglect the non-corotational drift motion. 

In the polar cap model of Ruderman & Sutherland (1975), the subpulse drift of the 
radio band is attributed to cE±/B drift in the polar cap region. The sub-pulse drifts have 
been observed from the old pulsars such as PSR B0943+10 ( P3/P1 ~ 37.35), B0826-34 ( 
~ 14), B2303+30 (~ 23), B0031-07 (> 34) (Gil, Melikidze & Geppert 2003). The pulsars 
quoted are old (r ~ lO^yr) and slow (P ~ Is) pulsars with spin-down energy many orders 
of magnitude less than known 7-ray pulsars (r ~ 10'^ ~ lO^yr). Even though the pulsar is 
old, the magnetic pair creation process in polar cap region must take place due to the high 
magnetic field {Bs ~ lO^^G). In the outer gap of such old pulsars, it is not clear whether the 
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high energy radiation and the photon-photon pair creation process take place efficiently. In 
following, therefore, we investigate the outer gap of the old pulsars. 

For simplicity, we assume that the available potential drop, (f)a = Bs^l^Rl/2c'^, on the 
stellar surface is used only for the acceleration of particles in the outer gap. This assumption 
is not valid for the young fast pulsars, but not bad for the old and slow pulsars because 
the gap covers the almost pulsar magnetosphere. The corresponding Lorentz factor will be 
estimated from mgC^Fa = e^^, to be ~ 2.5 x lO^P^^-^P^fg, where P is the stellar rotation 
period in units of Is, and P_i5 is the period derivative in units of 10~^^ s • s~^. The pair 
creation condition is E^Ex ~ 2(meC^)^, where Ex is the energy of the soft photons, and 
E^ is the curvature photon energy, which is E^ ~ O.lP'^'^PlfgGeV, where we adopt the 
light radius as the curvature radius of the magnetic field lines. From this condition, the 
threshold energy of soft photons in terms of P and P_i5 becomes Ex ~ SP^-^PTig^keV. 
For any neutron star cooling models and the polar-cap heating models, such high thermal 
temperature for old slow (P ~ Is, P_i5 ~ 1) pulsars is not predicted (Zhang et al. 2004, and 
references therein). For checking, we applied the present numerical code for an older pulsar 
PSR 0943+10 with the X-rays field of an unexpectedly temperature, which wc adopted for 
Vela pulsar [Ex ~ O.lSkev). We found that no pairs are produced in the gap. On these 
ground, we conclude that the pair creation process does not take place in the outer gap of 
the old pulsars quoted above. 

Such silent gap does not produce current carrying particles by itself. However if the 
external current comes into the gap through the boundaries, the electrons (or positrons) 
are accelerated in the gap, and then radiate curvature photons (for example, the energy of 
the 7-rays is about E^ ~ lOOMeV with P ~ Is and P_i5 ~ 1). If the external current has 
the Goldreich- Julian value Bg^l'^Rl/c, then the luminosity is same order of the spin down 
energy, namely Lj ~ Iqj x 0^ ~ PqD^P^/c^. For example, the predicted luminosity for 
PSR B0826-34, the observed distance {d ~ 0.54kpc, Lyne & Smith 1998) of which is the 
nearest in the pulsars quoted above, is ~ 6.2 • 10^° erg/s. However, since the corresponding 
flux on Earth, F ~ Lj/AncP, is ~ 2 • 10"^^ erg/cm^s, which is below sensitively of and 
unobservable even with the new 7-ray telescope GLAST. On the other hand, the middle 
age pulsars (r ~ lO^yr) will have the gap in which both of the pair creation and the 7-ray 
radiation processes take place, and are interesting candidates for GLAST (Torres & Nuza 
2003). 
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Figure 8. The non-corotational drift motion around the rotation axis on the outer boundaries. The solid, dotted, dashed, 
dashed-dotted correspond to the Case 1, 2, 3, and 4, respectively. The abscissa is the drift velocity in units of light velocity. 
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